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An equation proposed by Levy, Perdew and Salmi in 1984 [PRA 30, 2745 (1984)] is an orbital- 
free formulation of density functional theory. However, this equation describes a bosonic system. 
Here, we analyze on a very fundamental level, how this equation could be extended to yield a 
formulation for a general fermionic distribution of charge and spin. This analysis starts at the level 
of single electrons and with the question, how spin actually comes into a charge distribution in a 
non-relativistic model. To this end we present a space-time model of extended electrons, which 
is formulated in terms of geometric algebra. Wave properties of the electron are referred to mass 
density oscillations. We provide a comprehensive and non-statistical interpretation of wavefunctions, 
referring them to mass density components and internal field components. It is shown that these 
wavefunctions comply with the Schrodinger equation, for the free electron as well as for the electron 
in electrostatic and vector potentials. Spin-properties of the electron are referred to intrinsic field 
components and it is established that a measurement of spin in an external field yields exactly two 
possible results. However, it is also established that the spin of free electrons is isotropic, and that 
spin-dynamics of single electrons can be described by a modified Landau-Lifshitz equation. The 
model agrees with the results of standard theory concerning the hydrogen atom. Finally, we analyze 
many-electron systems and derive a set of coupled equations suitable to characterize the system 
without any reference to single electron states. The model is expected to have the greatest impact 
in condensed matter theory, where it allows to describe an iV-electron system by a many-electron 
wavefunction ^ of four, instead of 3N variables. The many-body aspect of a system is in this case 
encoded in a bivector potential. 



I. INTRODUCTION 
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Electronic structure simulations today underpin many 
models constructed to account for experimental data. 
Most of these simulations are based on an implemen- 
tation of density functional theory (DFT). DFT itself is 
theoretically justified by the Hohenberg-Kohn theorems 
[1], even though a model based on the electron density 
was suggested much earlier by Thomas and Fermi [2, 3]. 
A combination of Thomas-Fermi model and Hohenberg- 
Kohn DFT leads to the following description of a many- 
electron system [4], 
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-v H [p](r)+v ne [p}(r)+v xc [p](r) = p. (1) 



Here, T[p] is the kinetic energy functional, depending on 
the density of electron charge p, vh the Hartree potential 
or the electron-electron repulsion, v ne the electron-nuclei 
attraction, and v xc the exchange correlation potential, p 
is the Lagrange multiplier due to the condition of charge 
conservation, or the chemical potential. 

While this equation is generally valid, a transferable 
and fully general kinetic energy functional based on the 
density has so far remained elusive. It is known that the 
functional for the homogeneous electron gas is described 
exactly by the Thomas-Fermi kinetic energy functional 
(TF), or [2, 3, 5] 
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In present implementations of orbital-free DFT one ei- 
ther interpolates between the TF and vW functionals for 
more general density distributions [5, 7], or one bases the 
whole calculation or parts of it on the Kohn-Sham (KS) 
method of DFT [4, 8]. The reason, one tries to avoid the 
KS method as much as possible is that the computational 
effort in this case scales with the cube of N, the number 
of electrons, which limits the number of atoms which can 
be included in the simulation to a few thousand. How- 
ever, most of the technologically interesting materials e.g. 
in the semiconductor industry are very low doped; one 
dopant per millions of atoms of the host matrix. In this 
case current methods are unable to reliably simulate the 
material. 

In a paper in 1984 Levy, Perdew and Sahni [9] showed 
that one should be able to describe the general density 
of a system by: 
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For hydrogen atoms, it is described exactly by the von 



Here, v ex t is the external ionic potential, v e ff is the 
electron-electron repulsion combined with the exchange- 
correlation potential, and p is again the chemical poten- 
tial, which can be interpreted as a generalized and con- 
stant energy density. The problem, which was left un- 
solved, is that the equation seems to describe a bosonic 
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system. Subsequently, Norman March and others ex- 
tended the equation to account for the fermionic nature 
of electrons. They introduced a Pauli potential, which 
describes Pauli repulsion between individual electrons in 
a system [10-12]. 

This was the starting point of our analysis of orbital 
free DFT some years ago. The appeal of the equation 
is that it has the same general form as the Schrodinger 
equation and that it only includes the density, thus is an 
orbital-free formulation of DFT. But it is necessary to 
extend the equation and to generalize it, so that it can 
be used for arbitrary fermionic systems. This, as will 
be shown, can be accomplished in a consistent manner. 
While previous approaches focussed on the extension of 
the effective potential [10], within the present framework 
we extend the square root of the density to incorporate 
electron spin. The analysis necessary to accomplish this 
task is very fundamental, since it has to start with the 
question, how spin actually comes into a density distri- 
bution in the non-relativistic limit. Such a fundamental 
analysis has to include not only the structure of elec- 
trons, but also its interactions with electromagnetic fields 
to establish, that the model is fully consistent with the 
standard results. The structure of the paper is as follows: 

In Section II we review existing electron models in view 
of their suitability for the tasks at hand and give four 
postulates, which we think capture the requirements for 
a comprehensive electron model. We also introduce the 
wavefunction ip of single electrons as a multivector of 
three dimensional geometric algebra. 

In Section III we give a short overview over geomet- 
ric algebra necessary for the presentation. As condensed 
matter theorists will most likely not be acquainted with 
geometric algebra, we think this section is necessary for 
an understanding of the concepts. 

In Section IV we present a motivation why potential 
components of electrons in motion are necessary to attain 
wave properties of moving electrons. These electromag- 
netic energy components differ from the classical compo- 
nents previously introduced by Lorentz and Abraham. 

In Section V we present a consistent and local model 
of a single-electron's density and spin components in ge- 
ometric algebra as well as a relation between these com- 
ponents and the wavefunction of the electron. 

In Section VI we determine the interactions of the elec- 
tron with external static fields and show how the density 
and spin components change upon interaction. 

In Section VII we analyze the interaction of electrons 
and photons and show that the effect on the electron's 
properties can be mimicked by a calculation including 
only momenta and energies of electrons and photons. 

In Section VIII we show that the model is consistent 
with the standard model of hydrogen. However, we also 
present an extension of the standard model in the non- 
relativistic regime and including the spin of the electron. 

In Section IX we present the model of hydrogen 
molecules and extend Eq. (4) to an arbitrary density 
distribution and including the spin state of electrons. 



Finally, in Section X we discuss the presented concepts 
in view of fundamental issues in quantum mechanics and 
the consistency of the model with existing experiments, 
in particular the experiments on fundamental quantum 
mechanics. We also present an outlook on future work 
to be undertaken to check on the method. 

The guiding principle of this work is the 'radical con- 
servative - ism' of John Archibald Wheeler [13]: Insist 
on adhering to well-established physical laws (be conser- 
vative) but follow those laws into their most extreme do- 
mains (be radical), where unexpected insights into na- 
ture might be found. We do not attach any ontological 
claim to the model suggested. It is adopted purely for 
its usefulness, as it greatly simplifies the description of 
many-electron systems, in fact allowing for a coherent 
theoretical model from the level of single electrons to, in 
principle, macroscopic systems. 



II. SINGLE ELECTRON MODELS 

Generally, it has to be conceded that all attempts 
to ascribe a reality to the electron, which goes beyond 
a point particle with an intrinsic momentum and thus 
magnetic moment, have failed. The most convincing ar- 
guments against an internal structure and for a point- 
like electron come from single-electron diffraction and 
high-energy scattering experiments. In the first case, 
the point-like impacts of the single particles, gradually 
building up a diffraction pattern, are one of the great 
paradigms in quantum physics, which establish the exis- 
tence of a guiding principle behind the statistical events 
[14]. Similarly, high-energy scattering experiments are 
routinely used to check the predictions of quantum field 
theory. So far, the agreement between experiment and 
theory is better than one part in a billion [15]. 

There have been attempts, most notably by Louis de 
Broglie, David Bohm, and David Hcstenes [16-21], to 
ascribe the observed duality of the electron, its wave- 
features and its discrete charge and mass without mani- 
fest volume, to a field-like construct. The pilot-wave the- 
ory, originally due to Louis de Broglie, has even gained re- 
newed interest in recent years, see for example the article 
by Antony Valentini in a recent Physics World issue [17]. 
However, it is not clear, how these models can be recon- 
ciled with the fundamental fields of interaction known in 
physics. These interactions are all subject to relativity 
and thus do not allow for action-at-a-distance. Ascrib- 
ing a physical reality to the wavefunction of an electron 
itself, while still retaining the mechanical properties of a 
point particle, inevitably seems to lead to a pilot-wave or 
Bohm-type theory, where the potentials themselves be- 
come non-local [18]. It thus only shifts the central prob- 
lem, contained in the following question: " What actually 
is an electron, and which physical property describes its 
spin?" to the question: "What actually is a quantum 
potential and in which way does it relate to fundamental 
interactions?" In addition, it remains unclear, in these 
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models, how spin enters the description of electrons, and 
how it interacts with magnetic fields. As this is the fun- 
damental problem to be solved in the present context, 
one needs to go beyond the existing frameworks. 

If one takes current experimental results at the atomic 
scale seriously, then one arrives most likely at the con- 
clusion of Don Eiglcr [22] : " I don't believe in this wave- 
particle duality ... I think its mostly just the left-over 
baggage of having started off understanding the world in 
terms of particles and then being forced, because of the 
quantum revolution, to think of the world in terms of 
waves. Don't even think about them as particles. Elec- 
trons are waves. And if you think of them in terms of 
waves, you will always end up with the right answer. 
Always." A similar point of view, based on a mathe- 
matical analysis of quantum mechanical concepts, has 
recently been put forward by Hrovje Nikolic [23], call- 
ing the wave-particle duality a " myth" . In addition, our 
most successful theory to describe solids, density func- 
tional theory (DFT), even though it is based on second 
quantization, still ascribes a reality only to the density 
of electron charge and its field of interaction [1]. If this 
is the case, then a model of electrons could also be con- 
structed, not from mechanics and an elementary mass or 
charge, but from a wave-like structure which interacts in 
the same way electrons do. 

At this point one might observe that it is simply im- 
possible to account for the electron's spin within a clas- 
sical model, i. e. a model where spin is represented by 
a vector in real-space. However, as Hestenes, and Doran 
and Lasenby have shown [24, 25], the algebra of electron 
spin, described by the Pauli matrices, is also obtained by 
a vector model within Geometric Algebra (maybe better 
known as Clifford Algebra), where it results from the gen- 
eral properties of rotations in three-dimensional space. 
Along the same lines of thought one also arrives at the 
result that the 7-matrices in the Dirac equation of the 
electron can be seen as expressions of the geometric alge- 
bra of four-dimensional space-time. From these facts it 
is possible to conclude that spin can actually be seen as 
a "classical" geometric property of the spacetime repre- 
sentation of the electron, which can be described by the 
even elements of a multivector in geometric algebra [25]. 
If we choose to represent the wavefunction in a manner 
which bears on these relations, then it can be written in 
the following form: 

ip = a + ej A Gjfiij a, ,% e R (5) 

Here, the symbol A defines the "wedge-product" in ge- 
ometric algebra between the vectors of the three- 
dimensional frame, which yields a " directcd-planc" in 
space. Incidentally, in a frame of three dimensions this 
wedge product can be interpreted as the conventional 
cross-product x between two vectors times a "pseu- 
doscalar" , which in turn can be seen as the imaginary 
unit, or i. In three dimensional space the imaginary com- 
ponent can be interpreted in three different ways : (i) It 



is either a " dirccted-plane" , which means it is due to the 
directions of two independent vectors, or (ii) it is a vec- 
tor, since every directed-plane is dual to a vector, or (iii) 
it represents an angle and an axis of rotation [24, 25]. 
We choose the first interpretation in the following, be- 
cause it allows for a straightforward interpretation of the 
observed oscillations of electron waves as charge density 
oscillations. 

We shall start from some ideas developed more than 
ten years ago [26] , and show that they can be extended to 
yield a model of electrons which is free of contradictions, 
in accordance with most fundamental results obtained in 
quantum mechanics, shows that a "wavefunction" of the 
electron exists which must comply with the Schrodinger 
equation, and which allows a representation of electron 
properties within geometric algebra. The clear difference 
to previous models, e.g. the ones by Louis de Broglic [16], 
David Hestenes [20] or Jaime Keller [27] is that (i) the re- 
lation between wavefunctions and physical variables like 
mass density or field amplitudes is made explicit, and 
(ii) that it is shown in a local and time dependent pic- 
ture how external applied fields directly translate into a 
change of the wave-properties of the electron. 

The clear disadvantage of the model is that it is un- 
suitable to describe single-electron diffraction [14]. In 
this respect, the interpretation of the electron as a point 
particle and the interpretation of the wavefunction as a 
probability measure seems indispensable. 

However, the model does have conceptual advantages 
in that it allows for a direct and locally defined relation- 
ship between physical properties and (single-electron) 
wavefunctions. This gain, we think, makes it worthwhile 
to interpret electrons in this way. Moreover, it allows 
for a very efficient formulation of many-electron prob- 
lems. There, it leads to the introduction of a conceptu- 
ally new bivector potential. This potential is necessary to 
account for spin-properties in a correlated system. The 
description also remains remarkably simple and is thus 
potentially very useful in the development of numerical 
methods for the simulation of solid state and molecular 
systems. 

The model of electrons rests on four distinct postulates: 

1. The wave properties of electrons are a real physical 
property of electrons in motion. 

2. Electrons in motion possess intrinsic electromag- 
netic potentials which are vector-like. 

3. The magnetic moment of electrons is a consequence 
of the orientation of these electromagnetic poten- 
tials. 

4. In equilibrium the energy density throughout the 
space occupied by a single electron is a constant. 

It is a discerning feature of postulates that they cannot 
be derived from phenomena. The same applies to these 
postulates. However, their logical consequences have to 



4 



be in line with experiments or a generally applicable the- 
oretical framework. Moreover, they should make these 
consequences more transparent, i.e. they should lead to 
a gain in understanding. Finally they should enhance our 
ability to predict experimental results, i.e. they should 
lead to an extension of theoretical methods. We think, all 
this applies to our postulates, as we shall show presently. 



III. A SHORT INTRODUCTION TO 
GEOMETRIC ALGEBRA 

Most solid state physicists will not be familiar with ge- 
ometric algebra. This short introduction is modeled on 
an introductory paper by Gull, Lasenby and Doran [28] , 
which in our view provides the easiest introduction to the 
topic. In a Gibbs vector algebra, built on Euclidean ge- 
ometry, we have two separate products between vectors, 
an inner product a • b, which is a scalar: 

a • b = a e R, (6) 

and an outer product a x b, which is a vector: 

axb = cei? 3 . (7) 

In three dimensions the vector product is equal to the 
wedge product times the imaginary unit i, so that a x 
b = — i& A b. In geometric algebra the two products are 
combined in a Clifford product or the geometric product 
between vectors. 

ab = a b + a A b (8) 

The geometric product contains thus two parts: a scalar, 
which is symmetric, i.e., it does not change its sign upon 
a change of the order of a and b, and a bivector, which 
is antisymmetric: 

a ■ b = — (ab + ba) = b ■ a 

aAb = — (ab — ba) = — b A a (9) 

The product of two parallel vectors is a scalar, while the 
product of two orthogonal vectors is a bivector, denoted 
by the wedge product. The general product usually con- 
tains both, a scalar and a bivector. The wedge product 
can be seen as a plane, more specifically a plane, the ori- 
entation of which is given by the order of the two vectors 
in the wedge product. In geometric algebra, this is called 
a directed plane. 

It is important to realize that the geometric product, 
which adds a scalar and a bivector, will have properties 
of both, the scalar and the bivector. This is exactly like 
the relation between the real part and the imaginary part 
in a complex number. In two dimensions it can be shown 
that the wedge product carries the imaginary unit. The 
proof is simple. We assume two framevectors, ei and e 2 , 



which arc orthogonal. Then we get for the square of the 
wedge product: 

(ei A e 2 ) 2 = eie 2 eie 2 = -eieie 2 e 2 = -1 (10) 

This means that the wedge product carries an additional 
unit, the imaginary unit. A multivector, i.e., a sum of a 
scalar and a bivector in two dimensions can therefore be 
mapped onto the complex number space: 

z = x + yeie 2 = x + iy (11) 

At a deeper level this property is related to the fact that 
a wedge product is actually a rotator, or an element of 
geometric algebra which rotates vectors. This can be seen 
by acting with a wedge product on the frame vectors: 

(eie 2 )ei = — e 2 eiei = -e 2 

(eie 2 )e 2 = eie 2 e 2 = ei (12) 

Rotations are in most undergraduate physics courses de- 
scribed by matrices. The possibility to describe rotations 
by geometric products suggests a deeper link between ge- 
ometric products and matrices. In three dimensions a 
general multivector is described by a scalar, a vector, a 
bivector and a trivector. The geometric product between 
frame vectors has the following algebra: 

% = eie 2 e 3 
eie 2 = ie 3 

e 2 e 3 = ie 1 (13) 
e 3 ei = ie 2 
e l e i = 1 

Here, i is actually a pseudoscalar, or the highest geomet- 
ric product (the product of three independent vectors) 
in three dimensional space. The algebra of three dimen- 
sional space is equivalent to the algebra of the Pauli ma- 
trices, given by [28]: 

e i e j = % + itijkVk 

<r t dj = ISij + ieijkfrk (14) 

The Pauli matrices in this case can be seen as ma- 
trix representations of the 3-dimcnsional geometric al- 
gebra. A similar relation between framevectors of a 4- 
dimcnsional geometric algebra and Dirac's 7-matrices in- 
dicates that these matrices also can be seen as an expres- 
sion of 4-dimensional spacetime [25]. 

From a general point of view it thus turns out that 
imaginary components are quite frequently an indication 
of a geometric basis behind our formulations in a given 
vector space. It will be seen that this allows for a natu- 
ral distinction between density related and spin related 
components in the electron wavefunction. 

IV. WAVE PROPERTIES AND ELECTRONS 

In this section we review the original motivation to 
postulate the existence of field-like energy components 
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for an electron in motion. Wc show that if one assumes 
that wave properties of electrons are similar to the wave 
properties of electromagnetic fields, then the total energy 
of an electron will be different from its kinetic energy. As 
stated in previous publications, this is the only possibility 
to account for wave properties in a physical picture [26]. 
The approach in quantum mechanics, where one assumes 
from the outset that the only energy component of an 
electron must be its kinetic energy, leads invariably to the 
conclusion that the wave properties of the electron cannot 
be real which requires then to interpret them as related 
to the probability density [29] . This has been pointed out 
some time ago. However, here we start from the picture 
of waves developed in electrodynamics, and shall extend 
this picture to encompass also electrons in motion. In 
electrodynamics, the relation between the wavelength A 
of an electromagnetic wave and its frequency in vacuum 
v is given by the dispersion relation: 



\v = 



(15) 



Here, c is the velocity of a photon in vacuum. For photons 
we also must account for Planck's relation between the 
energy E and the frequency v of a particle: 



E = hv 



(16) 



Here, h is the Planck constant, or 6.626 x 10~ 34 Js. For 
electrons, we start with the de Broglie relation between 
wavelength and momentum of an electron. This rela- 
tion was first verified by electron scattering experiments 
of Davisson and Germer in the 1920s [30] and is today 
checked routinely in many labs around the world per- 
forming low energy electron diffraction experiments. It 
states that: 



A 



mv e i 



(17) 



m in this case is the mass of the electron, or 9.1 x 10~ 31 kg. 
To see, what happens, if we assume that the Planck and 
de Broglie relations are both valid for an electron, we 
now take the velocity c in Eq. (15) to be the velocity of 
the electron v e i, and combine it with Eqs. (16) and (17). 
Then we find that the energy of the electron should be 
twice its kinetic energy, since: 



A 



Vel 

E/h 



mv e i 



E 



mv 2 el = 2 x E k i 



(18) 



If we assume that an electron has wave properties which 
are similar to the wave properties of a photon, then we 
immediately arrive at the result that its energy is not 
its kinetic energy alone. It is, of course, not very plausi- 
ble that its energy is double its kinetic energy, since this 
would require in every energy balance of electron accel- 
eration or deceleration that a factor of two should mys- 
teriously show up in the balance sheet. This has never 



been observed, so the conclusion is fairly safe that this 
cannot be strictly correct. However, it illustrates a gen- 
eral point which will be the basis of our electron model: 
wave properties related to physical properties of single 
electrons (not a statistical manifold of many electrons) 
are only possible if the electron possesses more than just 
its kinetic energy. But it is also clear that the energy of 
the electron must be on average its kinetic energy, be- 
cause any other assertion requires too many additional 
assumptions to be reconciled with the energy principle. 



V. OSCILLATING ELECTRON DENSITIES 
AND POTENTIALS 

Concerning the electron, one could ask what property 
of an electron could actually be responsible for its wave 
features, if these are taken to be real. The obvious answer 
to this question, and one which will not seem strange to 
condensed matter theorists, is electron charge (or mass) 
density. In condensed matter these charge density waves 
are actually routinely observed, in particular in atomic 
scale experiments on metal surfaces. 

It is therefore also quite natural to assume that the 
wave properties will be related to some form of density 
oscillation. Here, we assume that the number density 
varies in the interval [0,1], where 1 indicates the density 
maximum. In the simplest case, that of a free electron 
traveling at a constant velocity v e i in z-direction, these 
oscillations are described by the plane wave: 



( +\ P° 



/4tt 

I — z — Airvt 



(19) 



The amplitude po will be subject to constraints, for ex- 
ample the condition that the density integrated over a 
certain volume equals unity. This condition, however, 
need not concern us at this point and it will in fact be 
shown that all fundamental results necessary can be de- 
rived without this normalization condition. It is impor- 
tant, though, to check that the ansatz is compatible with 
the limit of inertia. For vanishing velocity we get, using 
Eqs. (16) and (17): 



lim p 

U e( ->0 



(20) 



which is compatible with our basic assumption that the 
wave properties of electrons and their oscillating density 
are a direct consequence of their state of motion. An 
oscillating density of charge for a free electron is ruled 
out, if no external or internal potential energy compo- 
nents are present, since it violates the energy conser- 
vation principle. We thus have to introduce potentials 
to make up for the periodic variations. As the simplest 
case, these could be thought of as oscillating E or B 
fields. The reason for choosing vector fields rather than 
scalar ones lies in the one known additional property of 
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an electron: electrons possess an intrinsic magnetic mo- 
ment or a spin. Since such a property is incompatible 
with a scalar distribution of held properties, it has to be 
related in the simplest case at least to some kind of vec- 
tor field. Whether this is sufficient, has to be established 
subsequently by an analysis of interactions of electrons 
with external electro-magnetic fields. The two additional 
fields we introduce are an electric £ and a magnetic ~H 
field, which are thought transverse as in the case of pho- 
tons, but with double the wavelength and consequently 
half the frequency of the density oscillations. 



( — z — 2-Kvt - 
( — z — 2-xvt 



(21) 



The helicity in this way can either be positive - it com- 
plies with the standard right-hand behaviour - or neg- 
ative - it is the opposite. This feature introduces two 
possible groundstates of the free electron, which shall be 
later identified as its spin-up and spin-down state. The 
additional phase cf> has been added to account for energy 
conservation of electrons at the local level. 

Electromagnetic components to the electron's mass are 
not a new concept. They have also been proposed by 
Abraham and Lorentz in their " classical" models of elec- 
trons [31, 32]. Compared to these classical models, the 
present one is different in two aspects: (i) The shape of 
the electron, e.g. a point or a sphere, is not imposed from 
the outset. Such a shape would show up in high-energy 
scattering events in the scattering cross sections. While 
this is true for protons or, more generally, for atomic nu- 
clei, it has never been found for electrons. One reason 
could be that electrons actually do not possess a defined 
shape, but that their shape depends on the potential en- 
vironment, (ii) The electromagnetic components are not 
constant, but depend on the state of motion of the elec- 
tron. This is, as analyzed in previous sections, necessary 
to reconcile the assumption with the wave features of 
electrons. 

Energy conservation requires that the energy density 
at every single point of the electron is a constant. The 
energy density of the field components at a given point 
is: 



E 



field 



i 



(22) 



cos 



27r 

— — Z — ZTTl/t 

A 



The simplest solution to the constant energy density 
problem, which is also the simplest solution for the wave 
propagation of an electron is a phase shift of the fields 
by 7r/2, so that 



2 

E field 



(23) 



2 £ o£o + 2^°^ 



sin 



27r 

— Z — 2~KVt 



At this stage the introduction of fields with the given 
frequencies and wavelengths becomes clear, because we 
can eliminate periodic components of the total energy 
density with the help of the relation for the cosine at half 
angles: 

2cos 2 (x) = 1 + cos(2cc) (24) 
The kinetic energy density is then, in a first step: 



1 



PQVel 



1 + 2 cos 2 



2tt 



2nvt\ 



1 



(25) 



And with the following ansatz for the amplitudes of the 
fields: 



1 2 
^PoKl 



(26) 



We get for the total energy density: 



1 2 1 2 

Etot = ^PoVel + ^PoVel 



2tt 



2tt 



2 cos — z— 2-Kvt 



+ 2shr [ — z - 2-Kvt 



\pwli (27) 



In contrast to the statement in Eq. (18), the model leads 
to the result that the total energy is equal to the kinetic 
energy of the electron, as also assumed in quantum me- 
chanics. At a fundamental level, this is due to one feature 
of the ansatz: the wavelengths of field components and 
mass density components arc different. The electron's 
energy density is determined by its mass density at the 
groundstate, i.e. v e i — > 0, even though the electron wave 
is a physically real - i.e. with a physical property, the 
density, periodic in space and time - feature. The rela- 
tion between the frequency v and the wavelength A of 
this wave is the same as for a de Broglie wave, i.e., the 
group velocity is equal to: 



\mv\ 



/2h) 



dui 

dk d(mv e i/h) 



Vel 



(28) 



At this point the model captures at least four fundamen- 
tal properties of the electron in the conventional model. 
These properties are: 

• The wavelength of the electron wave is inverse pro- 
portional to its momentum (de Broglie). 

• The frequency of the electron wave is proportional 
to its kinetic energy (Planck). 

• The total energy of the electron is just the kinetic 
energy of its rest mass at every point of wave prop- 
agation (Energy conservation). 

• The total density of a free electron including the 
fields is a constant and equal to the inertial electron 
density (quantum mechanical plane wave). 
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So far all variables introduced, the density as well as 
the fields are thought to be physical quantities, i.e. they 
should be in principle measurable. But the ansatz also 
leads to the conclusion that a wavefunction, if defined 
on the basis of this model and in line with the required 
properties, cannot be a physical quantity. This is shown 
in the following section. 



A. Stability of free electrons 

Before, we have to consider an additional problem, 
which becomes imminent, as soon as the electron is con- 
sidered to be a structure with a finite extension in space 
and not a point-particle without structure: the problem 
of Coulomb repulsion. In DFT this problem arises due to 
the requirement that one electron does not interact with 
itself. Exchange and correlation potentials thus have to 
be corrected for this constraint [33] . The situation is dif- 
ferent in classical electrodynamics, where the electron, 
as a point particle, would carry infinite electrostatic en- 
crgy(see p. 751 of Ref. [34]). This has been the starting 
point of early attempts by Abraham and Lorentz to cir- 
cumvent the problem by stating that the electron's mass 
was purely electromagnetic [31, 32]. The electron radius, 
its classical radius, is then about 2.8 fm. However, this is 
still orders of magnitude larger than the largest possible 
radius inferred from high-energy scattering experiments. 
In relativistic quantum electrodynamics the problem of 
electron self-energy is removed by renormalization [35]. 

In the present context the problem is rather trivially 
solved (with, of course, the remaining problem of elec- 
tron self- interference, see the discussion), by defining a 
cohesive static potential (j) co h, depending on the electron 
density. Electron charge densities are typically one elec- 
tron per sphere with a Wigncr-Seitz radius r s of about 
two atomic units, as found in solids. We start from a 
spherical charge distribution of density po, with a radius 
of r s . Then the repulsion energy is given by the following 
integral (we use atomic units in the following derivation) 



W 



Po 



dr 



4nr 6 1 



-4nr 



o 

-jr sPo 



3 r 

4tt o \ 3 3 



The potential associated with this repulsion, the Hartree 
potential, is thus: 



V H = V H (r s ) = -r- 1 



(30) 



For electrons in metals the typical value for r s is about 
two Bohr radii. There is no reason to assume that the 
density for a free electron will be substantially different. 
Thus the additional cohesive potential for a free electron 
should be: 

V coh = « (r s = 2au) a -8.16eV (31) 

07% 



For a free electron the sum of both potentials must be 
zero. The repulsive and attractive fields are (Q is the 
volume of the Wigner-Seitz cell): 



Wcoh 

n 

<t>coh 



I 4n 3 

-poPcoh = -po -rrr s p 



3 

5r, 



3 

5r s 

^coh - 



= 



3 

5r s 



(32) 



Summarizing the results of this analysis we find that 
a model of free electrons, which refers wave properties 
to oscillating charge densities and electromagnetic po- 
tentials, can be in line with conventional experimental 
and theoretical models because the total energy of the 
electron is in both cases equal to the kinetic energy of 
the electron, and complies with the fundamental Planck 
and de Broglie relations. The repulsive Hartree potential 
within a single electron is corrected by a cohesive po- 
tential, which depends on the density of electron charge. 
The total electrostatic potential in this case vanishes, as 
required by all conventional theoretical models of single 
electrons. For the following derivations we use atomic 
units, where e = m = c = tl = 1. 



B. Why wavefunctions are not physical objects 

The requirements for a wavefunction ip or a spinor in 
the context of the Dirac theory are that the wavefunction 
can be written as a complex number (vector for spinors), 
that a duality operation ift — > exists, and that the 
product of ijj and ip^ is positive and equal to the num- 
ber density of the electron. In atomic units the number 
density is equal to the mass density and also, multiplied 
with the square of the velocity of light in vacuum, the 
energy density. In the following analysis this eases the 
notation considerably, while it would have to be com- 
plemented by suitable constants if performed in the SI 
system of units. Given the two separate energy contri- 
butions in our model, we can readily identify the real 
part of the wavefunction as the square root of the num- 
ber density, or p 1 / 2 . The imaginary part is more difficult 
to pin down. The requirement here is that it must be 
related to the magnetic properties of the electron, or its 
" spin" . The simplest ansatz is due to geometric algebra 
and the geometric product between vectors. Following 
our model we assume that the velocity vector of the elec- 
tron is parallel to e 3 . For ease of notation we introduce 
here a Poynting-like vector, where the electromagnetic 
energy flux is given by (in atomic units c = 1): 

S = £% = eie 2 £o'Ho sin 2 (~^~ z ~ 2tti4 j (33) 
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Then the energy density of field components is given by 
the following scalar: 



S 



ie 3 S 



2k \ 
S = So sin z ( — z — 2irvt\ 



So = £oHq 



(34) 
(35) 

(36) 



If we require that the real part of the wavefunction does 
not change under a transformation from positive to nega- 
tive helicity, and that the imaginary part is antisymmet- 
ric, then we can write the wavefunction in the following 
way: 



,1/2 



5 1 / 2 e 1 e 2 = p 1 / 2 + l 5 1 / 2 e 3 (37) 



The wavefunction in this case is a multivector composed 
of a scalar (the "real" component) and a bivector (the 
"imaginary" component). The duality operation changes 
the helicity of the electron; so exchanging the direction 
of £ and H we obtain: 

= p 1 ' 2 + S 1 ' 2 e 2 e 1 = p 1 ' 2 - iS^ea (38) 

For the product of ip and ijp we get consequently: 

^iP = iP^ = p + S (39) 

And if we set, as before: 

2 f 2lT 

p = po cos I — z — 2k vt 
2k 



S = So sin — z — 2nvt 



A 



So — Po 



(40) 



Then we obtain the result that the product tjj'ijj corre- 
sponds to the inertial number density of the electron: 



tp'ip = p + S = po = constant 



(41) 



which is the standard result in quantum mechanics and 
also in density functional theory. However, within the 
present framework wavefunctions also have a physical 
content, contrary to their role in the standard model 
where only their square has a meaning, that of a proba- 
bility density. One could pin down the difference by say- 
ing that even though the wavefunction is certainly not a 
physical object of the same reality as an electromagnetic 
field vector or a scalar potential, it contains physically 
real objects like fields and mass (or charge) densities - or 
rather the roots thereof. 

The multivector ip does not comply with the time- 
dependent Schrodinger equation, because the single time 
derivative of tp will yield the sum of a vector and an 
imaginary number (the "odd" elements of a multivector 
in three dimensions). In our view this bears on the fact 



that the model developed so far is not covariant. For ap- 
plications in DFT this is not relevant, as the LPS Eq. (4) 
is based on the time-independent Schrodinger equation. 
In order to compare the extended model of the electron 
to standard theory, we therefore define the Schrodinger 
wavefunction tps as a complex number, retaining the di- 
rection perpendicular to the field vectors e3 as a hidden 
variable: 



1/2 

Po 



I —z — 2-Kvt 



2tt „ 

| — Z — 2'KVt 



i 2tt 

p Q expz ( —z — 2/Kvt 



- .I/ 2 



(42) 



The kinetic energy operator or the time differential act- 
ing on ips then returns the energy eigenvalue of the free 
electron: 

, 1 4tt 2 1 



d I 
i— L Vs 



2 X- 



-tps 



V el^S 



E^js 



-i ujips = = Eips 



(43) 



dt 

It is interesting to note that the current density J 
as defined from the continuity equation and the time- 
dependent Schrodinger equation will be a function with- 
out local density variations. From the continuity equa 
tion 

i +v -* = dt 

we obtain the following result for J 



(J 

V • J = -Triples + V • J = 0, 



(44) 



J = 



2i 
2i 

2i 
Po 

PoVel 

poVei e 3 



{pf^VS 1 ' 2 - S 1 ' 2 ^^) 1 ' 2 
2ir „ \ a 

— Z- 2'KVt — 



(45) 



2tt 



2nvt\ 



2ixvt 



2n 



2tt 



2wut I 



e 3 



- 2-KVt ) 



2tt 



-Z — 2'KVt 



e 3 



which is equal to the momentum density of the inertial 
electron mass. It thus gives a consistent picture of the 
electron as an entity without any inner structure. 

To summarize the findings of this analysis we may say 
that even though this model electron is an entirely phys- 
ical object, all of whose properties are described be stan- 
dard physical quantities, we can only describe it as a 
scalar complex number, i.e. a wavefunction, if this wave- 
function itself does not have a direct physical meaning. 



VI. ELECTRONS IN STATIC EXTERNAL 
FIELDS 

While quantum mechanics, as shown in the previous 
section, describes an electron by a complex number, a 
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wavcfunction, which is not strictly speaking a physical 
object, it does so consistently by assigning a current den- 
sity to the propagation of the electron's inertial mass. 
Then, of course, an external electric or magnetic field 
cannot affect the density distribution of the electron in 
a physically transparent manner. Within the present 
context, where all components of an electron in motion 
are actually physical objects, we have to develop an un- 
derstanding, at a very fundamental level, what actually 
takes place if an electron is accelerated. Let us assume 
that this is due to an external potential </>. From a phys- 
ical point of view, four discrete processes are bound to 
happen: 

1. The electron velocity will change. 

2. The electron's density distribution will change. 

3. The electromagnetic field components will change. 

4. The intensity of the external field cf> will be dimin- 
ished due to energy transfer. 

The way to account for all four processes is captured 
in the following equation: 



dt 



(46) 



The sink in the scalar field 4> here accounts for the trans- 
fer of energy from the external field to the electron. It 
is easy to prove that the internal variations of the mass 
density or the field S do not change this equation. Wc 
know from the preceding derivations that 



S = Po — constant 



Then it follows that: 



S- 



P 







d , t dv e i 



(47) 



(48) 



Therefore Eq. (46) is valid without restrictions. Again, 
as for the current density, wc find that the density vari- 
ations are hidden in the combined effect on density and 
field contributions. Interestingly, this is only the case be- 
cause the time derivative of the density variations is of 
equal magnitude but opposite sign as the derivative of 
the field variations. 



A. The problem of frequency 



In particular, it is still unclear, how an external poten- 
tial actually influences the frequency of an electron wave. 
However, a closer look at the fundamental experimental 
results leading to the development of quantum mechan- 
ics reveals that the change of frequency with the inten- 
sity of electrostatic fields is already implicitly contained 
in Einstein's work on photoelectron emission. There, the 
kinetic energy of a metal electron and its frequency are 
related to the frequency of the incident photon and the 
intensity of the metal potential 4> m via: 



(49) 



Generalizing this result, it can be said that the frequency 
Lo e i of an electron in an electrostatic field, compared to 
the frequency of a free electron, shifts due to the 
existence of an electrostatic potential with: 



(50) 



Then the density components and the field components 
of the electron in an external potential (f> will comply 
with: 



i>s = p^' + iS 1 ' 2 



= P 1/2 
P 1/2 = /Tecs 

S 1 ' 2 



1/2 . 

Po sin 



2tt , , x ' 

— Z - (w - <p m )t 

2tt 



A 



z- (cjq - (j) m )t 



1/2 

Po expz 



2tt 



(w ~ 4>m)t 



(51) 



It should be clearly understood at this point that the 
frequency transferred to the electron is not the origi- 
nal frequency of the photon, but that this frequency is 
diminished, within the metal, due to the effect of the 
metal's electrostatic field. This diminished frequency is 
the frequency which the photon could actually transfer 
to the electron. The free electron, which is then mea- 
sured, has this diminished frequency and a wavelength, 
which corresponds to this frequency: a mechanism of en- 
forcing a wavelength, which we shall analyze in some de- 
tails further down when wc present a detailed model of 
photon-electron interactions. At this point it becomes 
clear that this wavcfunction, corrected for a change of 
frequency due to electrostatic potentials compared to 
the wavefunction of a free electron, complies with the 
general Schrodinger equation for electrons in a potential 
V = -6 m : 



At this point it is unclear, what the frequency of an 
electron wave after acceleration in an external poten- 
tial will be. We cannot use the conventional reasoning 
in quantum mechanics that the total energy is reflected 
by the frequency, while the wavelength is connected to 
the kinetic energy. This reasoning, which is one of the 
heuristic arguments to introduce the Schrodinger equa- 
tion, is only justified on the basis of classical mechanics. 



l dt^ S 
.di(js 



-"05 = UJotpS 



W - 4>rn)lpS 
1 



V 



(52) 



Again, it should be kept in mind that while the wave- 
function itself is not a physical object, but a multivector 
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in geometric algebra, its components, the real part and 
the bi vector or imaginary part, are related to real phys- 
ical properties. Here, we have assumed one direction of 
motion of the electron. It is, however, clear that the 
same general relation must apply to a general direction 
of motion. Eq. (52) is therefore generally valid. 

Considering, for example, an electron in a one- 
dimensional well V of length L, we find the usual result 
that k = 2n/L ■ n. In this case a superposition of left and 
right traveling waves corresponds to: 



-u> 2 



I 2 - le^S 1 ' 2 

(V> + + i~) = p 1/2 



(53) 



The density distribution measured, e.g. in Mike Crom- 
mie's quantum corral experiments [36], is then the real 
electron density and not the probability density. 



B. Electron in a magnetic field 



Following Doran and Lasenby [25] we could define the 
spin of an electron as: 



1 



-ipe 3 ip 



(56) 



Since tp = p 1 ' 2 + e 1 e 2 S' 1 / 2 and ^ = p 1 ' 2 - e^S 1 ' 2 the 
spin in this case is perpendicular to the field plane: 

4>e 3 ^ = pe 3 - (pS , ) 1/2 e 3 eie 2 + (pS , ) 1/2 eie 2 e 3 
— 5 f eie2e3eie2 

= (p + S)e 3 + (pS) 1 / 2 (eie 2 + e 2 ei) e 3 = p e 3 
1 



s = 2^° e3 



(57) 



The spin, defined in this manner, is a constant vector as- 
sociated with the direction perpendicular to the S plane. 
However, this is not compatible with the experimental 
results on Stern-Gerlach experiments [38], because an 
external magnetic field perpendicular to the electron's 
trajectory would lead to a vanishing magnetic moment. 
Calculating the equivalent expression for the direction 
e 2 , which is perpendicular to the direction of motion, we 
obtain: 



The acceleration of electrons in an external magnetic 
field is, conventionally, described by the Lorentz force 
equation: 



dv , _. 

Po ~dt = Po ( E + v X B ) 



(54) 



In geometric algebra, the same relation is expressed in a 
Lorentz covariant manner and using the Faraday multi- 
vector F (see p. 157 of Rcf. [25]): 



F = E + iB 



F ■ v 



(55) 



We have shown in the previous sections dealing with elec- 
tron acceleration that the electron's momentum change 
is described by the product of the inertial electron den- 
sity times the change of the electron velocity. This result 
is in accordance with the standard model, e.g. in density 
functional theory, where the electron does not possess 
an intrinsic and variable structure. Here, we only need 
to add that the same applies to the component of the 
Lorentz force which is perpendicular to the electron's ve- 
locity. A separate verification is therefore not required. 

The situation becomes more interesting, if we analyze 
the change of the electron's intrinsic structure in a con- 
strained trajectory. One can think, for example, of an 
electron contained in the potential well of a atomic nu- 
cleus, which has to follow the velocity vector of the nu- 
cleus. Then, the changes due to the external magnetic 
field must act on the electron's intrinsic properties as 
well. The important question, in the context of quantum 
mechanics, is how an external field will affect the elec- 
tron's "spin". However, we have not yet defined what 
the spin of an electron actually is in the present model. 



1/2 



ze 3 



S 1 ' 2 



e 2 \P 



,V2 



ze 3 



S 1 ' 2 



= nVa^^i/a _ „V2 p „,-^ 5 i/2 + ie3 s 1 l 2 e 2 p 1 ' 2 



= p 1 e 2 p ' - p ' e 2 ie 3 : 
ie 3 S 1/2 e 2 ie 3 S 1/2 

[p + S)e 2 + i{pS) 1/2 (e 3 e 2 - e 2 e 3 ) 



(58) 



Po 



[ —z — 2uit 



In the last lines we have used the anticommutation of 
frame vectors and the trigonometric relations for half 
angles. The directions ei and e 2 , perpendicular to the 
direction of electron motion, yield equivalent results: a 
constant vector in one direction and an oscillating vec- 
tor in the direction perpendicular to it. The spin vector, 
defined in the same manner as above, would then be: 



1 



2 



'Air 

e 2 + sin [ — z — 2u)t j ei 



(59) 



Here we slightly modify our electron model: we assume 
that the plane of electromagnetic fields is not perpendic- 
ular to the velocity vector, but that the velocity vector 
encloses an angle of 45° with its projection in the plane. 
The spin vector for this electron can then be described 
by: 



1 



:Po 



e 3 + e 2 1 . (Ait 

± = — ± — = sin — z — 2u>t ei 

y/2 V2 V A 



(60) 



The spin average over one period t =\jv will then be: 



s ± dt = ±-p 



1 e 3 + e 2 



V2 



(61) 
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It should be clear, though, that the previous derivations 
remain valid also for this modification of the electromag- 
netic plane. The reason for this is that the scalar and 
bivector components are separately accounted for and 
every cross product will vanish for a term ip^ijj. The ori- 
entation of the bivector plane is a free variable in the 
description, the energy conservation principle restricts 
only the phase of the wave, but not the orientation of 
the plane. According to postulate 3 of the introduction 
the magnetic properties of the electron should be due to 
the orientation of the fields, which in the previous equa- 
tion has been linked to the orientation of the spin vector. 
Then we may define a magnetic moment density of the 
electron by: 



fx = 7 ( S ) = 7 



p ±(e 3 + e 2 ) 
2 V2 



(62) 



In an external magnetic field, aligned for example in the 
e 2 direction, the magnetic moment density will interact 
with the field and lead to a variation of the static poten- 
tial according to: 



— — P ■ B e:I ; t — —ft ■ B n e 2 



(63) 



This static potential will affect the frequency of the elec- 
tron wave in the same way as an electrostatic potential. 
Writing, for efficiency, the wavefunction of the system in 
the magnetic field as: 



(64) 



where ipo encapsulated all other degrees of freedom of the 
electron wave and it is assumed that we can completely 
separate the magnetic degrees of freedom in the wave- 
function tp. Then from the time-dependent Schrodinger 
equation we get immediately: 



<9V 



UJ B 4> = ± 



-fB po 
2 V2 







± 



jB po 

2 yft 



(65) 



Within the framework of geometric algebra and using 
rotors instead of complex exponentials the wavefunction 
ip rotates around the direction of the external B field. 
The resulting effect is much the same: the interaction 
of the field-components with the external field leads to 
a rotation of the spin vector around the magnetic axis 
with two opposite vectors of magnetic moment [37] . In a 
Stern- Gerlach experiment, in this case, the ensuing tra- 
jectory will either be deflected along the positive or the 
negative y axis: exactly the result observed in the actual 
experiments [38]. 



C. Spin isotropy 

In standard theory, the spin of electrons is isotropic. 
This means, that the direction of the spin vector without 



an external magnetic field is undefined. It seems at first 
glance that the previous section, where the spin vector 
was oriented in a definite direction with respect to the 
vector of electron propagation, contradicts the model in 
standard theory. If, for example, the external magnetic 
field B is oriented along the ei axis, then the static po- 
tential (f>B will be zero. Clearly, this contradicts the ob- 
servations in Stern-Gerlach type experiments, where the 
split into two separate trajectories is independent of the 
orientation of the field. 

The aim of this section is to show, how the concept 
of isotropic spin fits into the presented model of elec- 
trons. From preceding sections we know that electrons 
with vanishing velocity do not possess field components 
to their energy density. If the plane of these field com- 
ponents is related to the direction of spin, then electrons 
of zero velocity do not possess spin. 

As the electron only possesses spin, if it is in motion, 
the question arises, how the spin vector relates to the 
plane of the electromagnetic energy components. The 
electric £ and magnetic H components are in the plane 
perpendicular to the velocity vector of the electron, the 
vector £ x % can therefore either be parallel (positive 
helicity) or anti-parallel (negative helicity) to the velocity 
vector. Then the spin vector of a free electron is described 
by the following relations: 



1 



1 



(66) 



In this case the spin-vector of a single electron is isotropic 
with respect to the plane perpendicular to the direction 
of motion. As this direction is singled out due to the mo- 
tion of the electron, the spin of a single electron cannot 
be completely isotropic, if spin is to be a consequence 
of electron motion. However, for a statistical manifold 
of n + spin-up (sg") and n~ spin-down (sq ) electrons, 
where n + = n~ the manifold will be fully isotropic in 
3-dimcnsional space. So far, the model is thus in accor- 
dance with experimental evidence. 

Considering the experimental evidence further, it is in- 
conceivable that the plane of the electromagnetic field 
components has a defined orientation, which is not per- 
pendicular to the vector of velocity, without any external 
B field, as this would lead not to two discrete experimen- 
tal results but, depending on the orientation of the field, 
to more or less continuous results. However, it is not in 
contradiction with experiments, if the orientation is well 
defined and not perpendicular to the velocity direction 
within an external B field. 

This analysis shifts the focus on the process consider- 
ably, because in this case the key problem is not to find 
an orientation of the spin vector, which yields two dis- 
crete results in every possible measurements in an exter- 
nal gradient of B (which, of course, will lead to the Pauli 
equation and Pauli matrices [39, 40]), but to determine, 
how an external field B changes the orientation of the 
field plane. If the external field rotates the field plane by 
an angle which is proportional to the angle between the 
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external field vector B and the vector of electron propa- 
gation v e ; , then the result is exactly the result obtained in 
the previous section: the magnetic moment then is either 
positive or negative with respect to the plane perpendic- 
ular to the external held and the possible measurements 
are either a deviation in positive or in negative direction, 
as observed in a Stern-Gerlach experiment. We have as- 
sumed, in the previous section, that this angle is half the 
angle between the field vector and the velocity vector. 
This choice, however, is not compulsory. As the exact 
angle will have to be determined by a careful analysis of 
experimental results, it is well beyond the scope of this 
paper, which only seeks to establish, that experimental 
results in static magnetic fields are in accordance with the 
proposed model of extended electrons. In fact, as shown 
in the following, this angle is contained in a constant, 
which describes the rotation of the field components by 
a modified Landau-Lifshitz equation. 
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FIG. 1: Rotation of the field vector of an electron described 
by a modified Landau-Lifshitz equation. 



D. Spin dynamics 

In condensed matter theory the dynamics of a spin sys- 
tem is described by the semi-empirical Landau-Lifshitz 
equation, which reads [41]: 



dM 
~dt 
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■MxH- 



(1 + cv 2 )A/ s 



Mx(MxH) (67) 



In this equation M is the magnetization vector, H the 
applied magnetic field, Ms the saturation magnetization, 
and a and 7 an empirical constant and the gyromagnetic 
ratio of the electron. As shown above, the field com- 
ponent S of the electron in motion will process around 
the magnetic vector if it is not parallel or antiparallcl to 
the velocity. This is described by the first part of the 
equation. The second part is commonly associated with 
damping. Here, we propose a formulation, which is sim- 
ilar to the second component of the equation, but will 
lead to two discrete and anti-parallel induced spin direc- 
tions. It includes the velocity vector of the electron to 
describe something akin to a torque acting on the field 
vector. The change of direction due to an external B 
field is described by: 



dS 
~dt 



const • S x v x 



~dt) 



(68) 



If the field is switched on in a finite interval, and to 
first order in the approximation, the induced field vec- 
tor within an external field is then described by: 



'induced 



= const • S x (v x B) 



(69) 



The equation has two discrete solutions (see Figure 1). 
Given now that the field vector precesses around the ex- 
ternal field we obtain exactly two solutions with opposite 
sign, which, in a field gradient along the x-direction will 
lead to two deflections in a Stern-Gerlach experiment. 



We conclude the presentation of a model of electron in- 
teractions with magnetic fields by pointing out that the 
key difference to standard models is that spin isotropy 
is taken to its logical extreme: the spin of an electron 
is only defined with respect to the velocity vector of the 
electron, but not with respect to the external magnetic 
field. External fields lead, as in standard condensed mat- 
ter physics, to a breaking of rotational symmetry. We 
may thus conclude that the model of electrons developed 
accounts for all known effects in static magnetic or elec- 
trostatic fields. 



E. Rotations and wavefunction symmetry 

The distinguishing property of fermions is that their 
wavefunction changes its sign upon a rotation by 2n. To 
see that this is also the case for our wavef unctions, we 
study a rotation of the spin vector s of the free elec- 
tron. In geometric algebra a rotation of a vector is ac- 
complished by a two-sided multiplication with a rotator 
R (Ref. [25], p. 274): 



s' = Rsltf 



(70) 



The one-sided rotator R for a rotation of 9 can be de- 
scribed as the exponent: 



B z = -1 



(71) 



R = R{9) =exp -5(9/2 
Since the spin vector is given by: 
e ± e ± 

the wavefunction must transform in the following way: 



(72) 



s' = R(6)ff> 
ip' = R(9)ip 



^R\9) 



(73) 
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That the wavefunction transforms via a one-sided and 
not a two- sided rotation is crucial for its symmetry prop- 
erties. Assume now that the angle of rotation 6 = 2vr. 
Then the spin vector will make a full rotation and re- 
vert to its original state. But this is not true for the 
wavefunction, since: 

ip = cxp —Birip = [cos(7r) — Bsin(ir)} ip = —ip (74) 

One could thus conclude that also the general symme- 
try of the wavefunction, which in many-electron systems 
is the basis of its construction from Slater determinants, 
comes from geometrical properties of the wavefunction 
as a sum of scalar and bivector terms. We shall use this 
property in our construction of the many-electron wave- 
function. 



F. Vector potential and wavelength 

An interesting variation of the same topic is the be- 
havior of electrons in an Aharonov-Bohm experiment. 
There, the magnetic field B is actually zero, while the 
vector potential A is non-zero. As measured orginally by 
Chambers [42] in 1960, and more recently by Osakabe et 
al. in 1986 [43], the vector potential affects the phase of 
the electron. The effect can be formalised as a change of 
the wavevector in an external field A with: 



k(A) = k + aA 



(75) 



Here, a is assumed to be a constant. That the potential 
must affect the wavevector or the local distribution of 
the electron density and the related potentials can be 
inferred from the description of a free electron after it has 
passed a region with non- vanishing potential A. In this 
case the wavefunction ips(z, t) corresponds to the original 
wavefunction ipso{z,t) augmented by an additional, and 
A dependent phase. For simplicity we assume that A is 
parallel to e3, then we get: 



fAs(A, z, t) = ipso{z, t) ■ exp ( ia / A(z)dz 

J z 



(76) 



Here, we have assumed that A is only non-zero from zo. 
Using the formulation of ipso for a free electron derived 
above, we may write for the wavefunction ips (zo is set 
to zero for convenience): 



ips(A,z,t) = expi 



2vr 



+ aA ) z — 2nvQt 



(77) 



Clearly the wavelength has been changed. If A is positive 
then the energy density in the region of the vector poten- 
tial has been increased. If this affects the electron, then 
it must lead to an increase of the wavevector. The con- 
stant a must therefore be positive. In atomic units and 
comparing with the experimental phase shifts one finds 
1. However, since the electron charge is negative we get 



for the wavefunction ip in an external vector potential the 
result: 



ips(A,z,t) = expi 



AO 



2-nvot 



(78) 



The density and field components of the electron in this 
case are described by: 



P 1/2 = P\ 12 cos 



S 1/2 = p 1 ' 2 



Pa sm 



V27T 


-*) 








V27T 






A Ao 





2-KVot 



(79) 



In an Aharonov-Bohm experiment with two different 
pathways characterized by opposite values of A within 
the interval zoj z ij the two partial waves at a point z will 
be: 



iP+(A,z,t) 



ip s (A,z,t) 



i ii 
p 1 cxp i 



1 12 

p 1 cxp I 



z- (zi - z )A ■ 
z+(zi — z )A ■ 



2-KVt 



2-KVt 



(80) 



A superposition of the two partial waves at z will conse- 
quently lead to an oscillating amplitude, described by: 

iPs = -(ip++ip-) 



i ii 
p 1 exp i 



2tt 

V 



z — 2-Kvt I cos(.zi — Zq)A 



(V /2 + iS 1/2 } cos(zi - z )A 



(81) 



This oscillation can then be measured, as in the case of 
electron density waves at surfaces, as an A and Z\ — zq 
dependent density variation at a detector screen. 

A generalisation of this result into a Schrodinger equa- 
tion which includes an external vector potential is not as 
straightforward as in the case of electrostatic fields. The 
reason is that A is not generally a vector, but a bivector, 
properly written as iA. (this follows from the fact that 
in a relativistic framework it has to be a complex field 
[25]). We start by writing the wavefunction in a more 
transparent form as: 



ips 



1 Ii 

p [cos ((k — iA) ■ x — cut) + i sin ((k — iA) ■ x — cut)] 

(82) 



Here, x is a three-dimensional position vector, and i is 
the pseudoscalar eie2e3. If we assume that the vector 
potential in the region of interest is constant, then a first 
derivative of this wavefunction leads to the following re- 
sult: 

Vi/>s = Po /2 {-(k-iA)sin((k-iA)-x-wt) 

+ (ik + A)cos((k-iA) -x-wt)} (83) 

k sin 4> + A cos <fi + i (k cos <p + A sin <p)] 



1/2 
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For brevity we have written </> for the argument of the 
periodic functions in the last line. It can clearly be seen 
that the inclusion of a general vector potential into the 
Schrodinger equation requires that it has a different na- 
ture than the terms related to the wavevector k. The 
gradient of ips in this case is a multivector composed of 
even (the bivector terms) and odd (the vector terms) el- 
ements. For the second derivative one obtains with the 
same procedure: 

V 2 Vs = Po 2 (-k 2 - A • k - k • A + A 2 ) Vs (84) 

To account for the changes in the local derivatives com- 
pared to the frequency in a Schrddinger-type equation 
one consequently has to set: 

^ = ^V-A) 2 <A S (85) 
Including also a scalar potential V , we obtain: 

'^-K*5r A )'*" +v *' (86) 

The equation is equal to the general non-rclativistic 
Schrodinger equation for electrons in the presence of elec- 
trostatic and vector potentials. The suggested extended 
model of electrons is thus fully compatible with the non- 
relativistic theoretical framework of quantum mechanics. 

VII. ELECTRONS IN DYNAMIC EXTERNAL 
FIELDS 

Based on these findings one may equally analyze the ef- 
fect of an external electromagnetic field on the velocity of 
electrons. The electromagnetic field shall be described by 

1/2 

a bivector iS e ( n e3, which propagates in e3 direction with 
velocity c, the velocity of light in a vacuum. The bivec- 
tor iSlm^z encaptures electromagnetic £ and % fields 
perpendicular to 63. Since there is no real component of 
the electromagnetic field vector, the field can also not act 
directly on the real part of the electron's wavefunction, 
or its density distribution. However, we found in the 
analysis of the electron's acceleration in an external field 
that the time derivative of the field components and the 
density components are equal and of opposite sign. One 
may infer from that finding that an external and oscillat- 
ing field S will lead to an oscillating field p via something 
akin to the Newtonian principle: for every action (here, 
a variation of field intensity due to an impinging electro- 
magnetic field), there is a reaction (here, a variation of 
mass density of the electron). To see whether such an as- 
sumption is justified we consider an electromagnetic field 
propagation in ^-direction with the field vector S em given 
by: 

S em (z,t) = 5 (i)sin 2 [^-z-2nut] (87) 



Given that the field propagates with relativistic velocity 
we assume, for this treatment, that the electron veloc- 
ity is very small and that in effect the propagating EM 
field interacts at a given point with the electron, which 
perceives the presence of the field more or less as a time- 
dependent variation. At a point zq, <po = {2n/X)zo this 
means: 

S em {z ,t) = S (t) sin 2 (fo - 2<irvt) (88) 
— AnvSv(t) sin (</>q — 2m>t) cos (</>o — 2itvt) 

If the effect of the electron's density is equal to the change 
of the external field 

dp e i{z ,t) = ~S em {z ,t), (89) 

then the density change at zq, which we assume to be 
initially constant and po, will be: 

dp e i{z 0l t) = po - S (z ,t) sin 2 (0 O - 2nvt) (90) 

According to the relation between density and field com- 
ponents this also leads to a change of the internal poten- 
tials S e i of the electron, described by: 

dSei(zo, t) = S (z 0l t) sin 2 (cf> - 2-Kvt) (91) 

The limit of energy adsorption from the electromagnetic 
field is reached, when the electron density has reached 
the limit: 

Pei(zo,tf) = Po - S (z ,tf) = (92) 

While this seems a coincidence at first view - the ampli- 
tude So could in principle be arbitrary and the electron 
adsorb less energy than required to reach this threshold 
on reflection it is quite understandable. After all, the 
external field S p h will have been created by the reverse 
effect, i. e. an electron being decelerated and emitting 
a photon in the process. Such an emission process is 
equivalent to an adsorption process with a negative time 
coordinate. From a physical point of view it is therefore 
symmetric. 

However, it has to be acknowledged that this partic- 
ular point deserves a more careful analysis, which has 
to be given in future. The second point, which has not 
been treated is the transfer of frequency, due to the exter- 
nal S field, to a change in wavelength, or a local change 
of electron density and field amplitudes after adsorption 
of energy. So far, this point has been omitted since we 
have analysed energy transfer processes at one particu- 
lar coordinate zo only. To determine the phase difference 
between different points at the moment, when the elec- 
tron acceleration is terminated (t = tf), we calculate the 
density at the point z + dz. From the equations derived 
so far it is clear that: 

7T 

(j>o-2'Kvtf = - S (z ,tf) = po (93) 
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Calculating now the density at zq + dz and t 
Pei(zo + dz,tf) = S (z ,tf) 



S {z + dz,tf) 



sm 



7T 2TT , 

2 + T dz 



So(z ,tf) — S (zq + dz,tf) cos 



Setting now: 



S (z + dz,t f ) w S {z ,t f ) = po 



we obtain for the amplitude of the density at the point 
zo + dz the result: 



tf we get: e. in an environment, where the electrostatic potential 
varies within the mass distribution of a single electron. In 
(94) particular, since we expect in such an environment that 
the frequency changes continuously. Within the present 
framework the Schrodinger equation, which describes the 
energy levels within a hydrogen atom, is due to a slightly 
different physical situation than treated previously. This 
situation is captured in postulate 4: In equilibrium, the 
energy density throughout the space occupied by a single 
electron is invariant. It can be argued that such a pos- 
tulate is necessary to account for the - experimentally 
observed - stability of hydrogen. Let us assume that the 
postulate is not valid. We may then write the equation 
in the following way: 



(95) 



p e i(z + dz,tf) = /9 sin 2 



dz 



(96) 



The local distribution of the density is then equal to the 
local distribution of a wave with velocity v = h/mX. A 
similar result is obtained for the field component of elec- 
tron propagation. Since both components show the same 
properties as an electron wave in motion with velocity v 
it is safe to conclude that the adsorption of the S field 
leads to an acceleration of the electron until it reaches a 
velocity which is corresponding to the frequency of the 
incident dynamic field. 

As both components comply with the standard form 
for the wavefunction, we may write the electron wave- 
function after acceleration as: 



^s(v) =p x > 2 {v)+iS 1 ' 2 {v) 



po exp i 



2tt 



(97) 

One may now analyze the problem of an electron in mo- 
tion undergoing accelerations due to impinging photons. 
Given that the wavefunction is linear, a subsequent ac- 
celeration will change it in exactly the same manner as 
the first one. This means, that one can actually analyze 
electron-photon interactions from the viewpoint of the 
conservation of energy and momenta. This is, how one 
commonly accounts for Compton scattering [44] . 



VIII. THE HYDROGEN ATOM 

It is straightforward to apply the Schrodinger equation 
of the electron, developed in previous sections, to the 
problem of a central potential and to restrict the problem 
to finding solutions of the equation: 



1 



1 



eips 



However, even though we have shown how the frequency 
of the electron changes in an external electrostatic field, 
we have done so only for one particular value of the elec- 
trostatic potential. We have not analysed, whether the 
electron can have a single energy value, described by 
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--V 2 



- ) 



e(r)^c 



(99) 



If a particular infinitesimal density component dp is 
found at a radius r, its energy value is described by e(r). 
Now if e(r) varies with the radius, for example becom- 
ing more negative with decreasing r, then the system can 
lower its energy if dp migrates to a position closer to the 
nucleus. In this case the general principle of energy mini- 
mization requires that the electron density collapses to an 
infinitesimal shell around the nucleus. Since this is not, 
what we observe, postulate 4 must be correct. And in 
this case the energy value e must be constant throughout 
the system. 

The standard solution for the non-relativistic hydrogen 
problem is the following (we give its standard textbook 
format): 



ipnlmir,^,^) = C„ im exp (-- — 



2r 
nag 

x P ; m (cosi9) exp(imip) 



r2Z+l 



( 2r 
\na 



(100) 



Here, C n i m is a normalization constant. Since the hy- 
drogen atom is neutral, the total electron charge within 
its shell must be equal to unity. a is the Bohr radius, 
or 5.292 x 10 _11 m. Writing the wavefunction in a more 
compact manner we may state: 



Ipnlmir,^,^) = U n l m (r,'&)e 



! Ili r ~ 



(101) 



Density components, as analyzed in great detail in the 
preceding sections, are reflected by the real components 
of electron wavefunctions, while the imaginary compo- 
nents relate to the electron's fields. This allows to de- 
compose the wavefunction of the hydrogen electron into: 



p l/2 = U n i m (r, f9) cos mip 
S l/2 = U n i m (r, tf)smmcp 



(102) 



Here, the important result is that in general the electron 
within a hydrogen atom also possesses density compo- 
nents and field components, related to its state of motion. 
Doran and Lasenby have given a fully relativistic treat- 
ment of the hydrogen problem within geometric algebra, 
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based on previous work of Arthur Eddington [45, 46], and 
found that the general solutions have to be described by 
multivcctors. This indicates the same feature: a presence 
of density components and field components (see p. 294 
ofRef. [25]). 

The field plane is perpendicular to the direction of e v . 
It should be kept in mind that this wavefunction de- 
scribes an electron without spin: the components ra^O 
relate to additional motion of the electron around its 
z-axis. Including spin into this description requires to 
modify all states, also the groundstate m = 0. As the 
groundstate solution of hydrogen is radially symmetric, 
and as the spin-component of the wavefunction is de- 
coupled from the local component and has exactly two 
solutions, we may account for spin by the component: 



+ 1 / 
X = ^/f ( 1±e,J Ae ^ 



(103) 



For the groundstate wavefunction m = the result will 
thus be: 



1 



^nlo^'ti'Viv) = -7=^nio(r,i?) (lie^Ae^) 



Writing the local part of the wavefunction as a multivec- 
tor leads to the following formulation: 

ipnim = U n i m (r, i9) (cos mtp + e r A e# s'mmtp) (104) 

The complete solution including spin then can be written 
as the following multivector: 



Vnlm — 



(105) 



V2 

x {cosmip + i [(e v =p e$) sin my =F e r cos mtp]} 



IX. THE HYDROGEN MOLECULE AND 
DENSITY FUNCTIONAL THEORY 

A hydrogen molecule or a helium atom poses the addi- 
tional problem to reconcile postulate 4 with the existence 
of one coherent electron density throughout the system. 
In principle, the problem can be solved using standard 
DFT based on the Kohn-Sham equations [8]. However, 
the Kohn-Sham formulation of DFT seems, from the out- 
set, somewhat cumbersome, since it requires the self- 
consistent calculation of every single-electron state in a 
system. The computational effort required as the system 
size increases then rapidly becomes prohibitive. This has 
led to extensive research into the possibility of a so-called 
orbital- free formulation of DFT [47] . The key problem in 
this line of research is the kinetic energy functional, which 
to date cannot be described in a completely transferable 
manner. Given the preceding presentation, it seems quite 
unlikely that a relation, based on electron density alone, 



could actually account for all aspects of a many-electron 
system. If magnetic properties play a role at all, then 
according to this model they must enter the theoretical 
description. This applies not only to magnetic material, 
but also to non-magnetic solids, as the minimum energy 
state of the system is a consequence of the existence of an 
equal number of spin-up and spin-down electrons. How- 
ever, the many-electron problem was encoded by Levy, 
Perdew and Sahni (LPS) in 1984 in an equation based 
only on the density [9] (see Eq. 4). 



A. Modified density equation including spin 

Within the present context, there are two problems 
with the equation: (i) It describes a bosonic system, since 
all properties, related to the spin of electrons, are absent 
from the formulation, (ii) It does not account for the 
difference in nature of field components (a bivector) and 
mass components (a scalar). It is hard to see, therefore, 
how this relation could be complete. However, it should 
be possible, given the analysis contained in this paper, 
to accurately describe a non-relativistic many-electron 
system in its groundstate by the following coupled equa- 
tions: 



Vext 



' J eff,0 



Pa 



(106) 



v e ff,o combines the Hartree potential and the cohe- 
sive potential defined in Section V, but does not con- 
tain exchange-correlation potentials. Exchange and cor- 
relation potentials encode the difference between non- 
interacting and interacting electrons. In this equation 
electron interactions beyond v e f f t o are part of the bivec- 
tor potential ie v Vi, since they are due to field- mediated 
interactions within the electron distribution. The direc- 
tion of the unit vector e„ and the intensity of the po- 
tential Vi will depend on the system under consideration 
and have to be determined self-consistcntly. The bivec- 
tor part of the many-electron wavefunction, the product 
of unit vector eg and electromagnetic intensity S 1 ^ 2 are 
equally system and location dependent and have to be 
determined by self-consistent iterations. The equation 
couples field terms to density terms also via a bivector 
potential. The resulting mass and field components of 
the wavefunction then give the conventional mass den- 
sity po- Even though both es and e v will vary from one 
point to the next, the set of coupled equations is still 
much simpler than the full many-body treatment of a 
many electron system in standard theory. 

The justification to write the wavefunction of a many- 
electron system in this way rests on three results of pre- 
vious sections: (i) The wavefunction of any single - and 
extended - electron complies with a Schrddinger equa- 
tion. A system composed of an aggregate of more than 
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one electron then should also comply with a Schrodingcr- 
like equation, (ii) The construction of a wavefunction 
from scalar and bivector components ensures that the 
47r-symmetry of the wavefunction will be retained, which 
is the key property of the wavefunctions of fcrmions. (iii) 
The interaction between field components and mass com- 
ponents implies bivector potentials in a crystal field of in- 
teractions. Without such a potential many-body effects 
could not be included in the description. 

One can decompose the equation into scalar and bivec- 
tor components to make transparent, in comparison with 
the LPS equation, what the additional components actu- 
ally are. For the scalar part we get: 
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+ V ext + V e fffi - M 



,V2 



Vl S 1/2 (107) 



The change is thus not simply a potential term of the 
same general form as the effective or external potential. 
Such a behavior could be encoded e.g. in a density de- 
pendent exchange correlation potential. However, here 
we have an additional term which will act as a source for 
the electron density distribution. This source term de- 
pends on the direction of the bivector potential as well as 
the field components of electron motion. It is thus clearly 
beyond current models in DFT. The bivector equation is 
equally instructive. It can be written: 



1. 



-V 2 + V ext + V e fffi - /U 



= e v A e 3 ViS 1/2 



ie v Vip 



ie s S 1 ' 2 

1/2 



(108) 



This is, in essence, a vector-equation with purely imagi- 
nary variables. It thus introduces directional effects into 
the density equation, and couples density components 
and field components. Also this equation is clearly be- 
yond current formulations in DFT. In imaginary com- 
ponents and replacing the wedge product by the cross 
product (a A b = ia x b), we get: 



1. 



-V 2 + V ext + V e fffi - M 



e„ x e s 



zuS 1 ' 2 



e v Vip 



e s S^ 2 

1/2 
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It can also be seen that either the field component S 
or the bivector potential Vi by itself will not alter the 
general equation. Only the field components of electron 
energy and the bivector potential of the solid together 
could give an accurate account of the charge distribution 
for electrons with spin. It is also clear, by comparing with 
the hydrogen problem, that the potential Wj is absent 
in the groundstate of the hydrogen electron: it is thus 
genuinely of many-body origin. 

In a slightly more compact form, which differentiates 
between vector and scalar components of the equation, 
we may write the following coupled set of equations, 
where the local dependency of components has been 



made explicit for clarity: 

A(r) - egS 1 ' 2 



II(r) 



~-V 2 + v (r)-p 
~V 2 + v (r)-p 



G v Vi 

Vext + Vefffi 



(110) 



p 1 / 2 ^) = II(r) • A(r) (111) 
A(r) + p 1/2 (r)n(r) = n(r) x A(r) 



In general, the formulation is still much simpler than 
present formulations of many-body theory. The simpli- 
fication comes from the fact that the model wavefunc- 
tion contains only four (the amplitudes p 1 / 2 and S 1 ^ 2 
and the vector of unit length es) independent variables, 
which can be mapped onto a grid in real space. It is 
thus much easier to evaluate in a general minimization 
problem even though it captures all physical properties 
of the system. Using, for example, a real-space method 
and a grid spacing of about 0.2 to 0.3 A, it should be- 
come computational routine to simulate systems of more 
than one million atoms: an increase of the current limit 
of DFT of a factor of about one thousand. In this case, 
mesoscopic systems will come into the range of ab-initio 
methods and biological systems should be tractable by 
theoretical means. 

The indispensable first step for such a theoretical 
framework will be to determine the bivector potential and 
the effective potential in a simple system, say a system 
of constant electron density. Only after such a poten- 
tial has actually been determined can the framework be 
applied to more complex systems with variable electron 
density. While it cannot be guaranteed that this ap- 
proach will lead to a reliable method for computing the 
physical properties of many electron systems, it seems 
justified to develop its consequences and to compare the 
results to standard methods. The gain in theoretical ef- 
ficiency, if successful, would be quite substantial. 



B. Many-body wavefunction 

The groundstate wavefunction of an A^-electron system 
can be written as: 



MB 



* (xi,x 2 ,x 3 , ...,xjv) 



(112) 



Here, the variables Xj combine local (r^) and spin (cr,) 
variables. The wavefunction ^mb thus contains at least 
3N variables, the position vectors of all N electrons, ri to 
i"tv ■ The main condition imposed upon the wavefunction 
due to the Pauli principle is that it must be antisymmet- 
ric: an exchange of two variable (or two electrons) and 
Xj (i and j) changes the sign of the wavefunction from 
plus to minus and vice versa. A general form of the wave- 
function complying with this condition is its construction 
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via Slater determinants [48]: 



= ^2 C K Det\^ 1 (x 1 ) 1 ^2{x 2 )^d,{x 3 ,), ...,ip N (x N )\ 

K 

(113) 

For a hydrogen molecule the many-body wavefunction of 
the system, including the spin functions a and /3, is given 
by [49]: 



^(ri,r 2 ,CTi,CT 2 ) = C 



^ioo(ri)o!(cri) -0ioo(r2)a((T2) 
■0ioo(i"i)/3(cri) ■0ioo(i"2)/3(o-2) 



The total (groundstate) wavefunction is therefore: 



(114) 



*(ri,r2,eri,cr 2 ) = ^100(^)^100^2) 

x \a{ai)p{a 2 ) - P{ai)a{a 2 )]{llh) 

Since a hydrogen molecule in its groundstate does not 
interact with a magnetic field, one may omit the spin- 
components of the wavefunction and write: 

*ff 2 (r) =Vioo(r-Ri)Vioo(r-R 2 ) = pj /3 (r) (116) 

where we have assumed that the nucleus of the first atom 
is at Ri and the nucleus of the second atom at R2. How- 
ever, including spin becomes necessary if one considers 
the gradual transition from the hydrogen molecule to two 
hydrogen atoms [50]. In this wavefunction we may ac- 
count for spin in a similar way as for the hydrogen atom, 
with the only difference that now we have two centers 
which serve as the origin of radial unit vectors. These 
unit vectors arc symbolized by e* and e 2 r , respectively. 



*ff 2 (r,a) 



-?/>ioo(r - Ri)Vioo(r ~ R2) 



[(1 



'1 



1 



1 
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In its triplet state the spin component of the two electron 
wavefunction is symmetric, while the local component 
becomes antisymmetric. In this case one has to construct 
the charge and spin densities from the Is and 2s orbitals, 
respectively. 



X. DISCUSSION 

We have on purpose avoided any reference to classi- 
cal mechanics in this work on electrons. The reason is 
that we think that mechanics has obscured the mean- 
ing of certain formulations in quantum mechanics. In 
that respect, we also think that the name quantum me- 
chanics carries a wrong message, and that it is better to 
talk about either wavedynamics or microdynamics, al- 
though this might be a secondary issue. The primary 
issue, we hope to have clarified in this paper, is that elec- 
trons are not mechanical objects. This, of course, is well 
known. However, so far no detailed account existed, how 
density and spin properties of single electrons might be 
distributed in space and time without any measurement. 



We think that the concept of an extended electron, intro- 
duced here, might be easier to reconcile with experimen- 
tal results at the atomic scale, with density functional 
theory, and with visualizations of physical processes. It 
is often the picture in the mind, which leads to new in- 
sights, and not the mathematical formalism. Whether 
this electron model describes the real electron or not, is 
a moot point. As demonstrated, it leads to the same 
numerical results as the standard model, it accounts for 
most experimental results obtained so far, it gives a con- 
sistent interpretation of the electron wavefunction, and 
it allows to picture atomic-scale processes in space and 
time. Even if it were not real, this electron seems quite 
a useful model. 

It is instructive to look at the presented material from 
the viewpoint of basic principles. Quite apart from the 
historical development of modern electron theory, this 
should allow us to appreciate the importance or unim- 
portance of specific findings. The first fundamental find- 
ing, and which made a wave theory of electrons initially 
necessary, is actually electron diffraction. That this was 
initially a hypothesis in a PhD thesis of Louis de Broglie, 
makes it all the more remarkable. In a modern context, 
one could say it is the observation of standing waves on 
a metal surface, which leaves no other conclusion but to 
assign wave properties to electrons. In this respect it is 
hard to overestimate the importance of the experiments 
coming out of Don Eigler's lab in the early 1990s. 

As a second fundamental finding one has to nominate 
photoelectric currents as interpreted by Einstein with the 
help of Planck's hypothesis. Here, it is found that the 
frequency of the electron waves can be altered by elec- 
trostatic fields, as existing within a metal crystal. This 
finding paves the way for the formulation of the initial 
Schrodingcr equation, which can be used to determine 
the emission frequencies of hydrogen atoms, admittedly 
one of the great problems of 19th century physics. How- 
ever, this finding alone does not allow for a treatment of 
interactions with electromagnetic fields. 

Here, we need the third fundamental finding, which 
we consider to be the Aharonov-Bohm effect. Only in 
this context can it be understood that the vector poten- 
tial changes the wavelength and thus has to be treated 
differently than, e.g. the changes in an electrostatic po- 
tential, where the frequency is affected. It is also quite 
interesting that this effect was a theoretical prediction 
based on the Schrodingcr equation and formulated long 
after the equation had been used to solve the hydrogen 
problem. The Aharonov-Bohm effect is as important as 
the photoelectric effect. Without it, the only way to in- 
clude vector potentials in the Schrodinger equation would 
be via a classical Hamiltonian H = [p — (e/c)A] 2 [34], 
and the correspondence principle. This, in turn, would 
require to treat the electron as a classical point-particle 
and lead to inconsistencies in the suggested theoretical 
model. 

The most intriguing part of the framework seems 
to concern the possibility of constructing wavefunctions 
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from density and field components and, by the reverse 
process, of decomposing wavefunctions into density and 
field components. This should add substantial physical 
insight to any theoretical result obtained e.g. within a 
many-body framework. In the final analysis we think 
that the role and the actual - i.e. physical - content of 
the wavefunctions have not been sufficiently clear. An 
example of this lack of clarity in our view is the famous 
Schrodinger-cat thought experiment [51]. There, a su- 
perposition of two states, a living cat and a dead cat, 
collapses into one final state upon measurement, i.e. the 
interaction of an electron with an external field. The 
important feature of the experiment lies in its statistical 
interpretation of the wavefunction. As we do not know, 
which state the electron is in, the cat could be said to 
be dead and alive at the same time - which seems quite 
impossible for a living organism. However, as shown in 
the treatment on spin and its interaction with a magnetic 
field, the electron wavefunction within the present frame- 
work will always indicate a well defined physical state, as 
described by time and space dependent density and field 
components. The statistics in this case do not enter via 
some unknown "quantum" property of electrons, but are 
due either to a variation of the electron's properties in 
space and time (which, for example, prevents us seing 
the density variations of a surface-state electron if it is 
not scattered on a defect, and makes us detect standing 
waves at a step edge or a single atom [36]), or to a lack 
of knowledge about the individual electron (which makes 
us detect two deflection spots in a Stern-Gcrlach exper- 
iment [38]). In this case, if the life of a cat is related to 
one particular electron state, then this cat is definitely 
dead or alive, but never both. If we have to describe 
it as a superposition of alive and dead, then we have to 
concede that we do not know exactly which electron we 
are talking about, or at which precise moment our mea- 
surement is taken. This is, in essence, the same lack of 
information encountered in statistical thermodynamics. 
Incidentally, this view is very close to the one advocated 
by Erwin Schrodingcr [51]. 

However, we have to concede that the model does not 
provide a comprehensive theory for the observed self- 
interference of single electrons [14]. Interference here is 
due to the spacetime variations of the electron's wave- 
function and its geometric characteristics upon scatter- 
ing. How a single electron could change its wave pattern 
by interacting with an environment which is well sepa- 
rated from any one of its possible trajectories, cannot 
be answered within this model. The only way, it seems, 
that such a behavior could be reconciled with the present 
model, is by invoking a detection loophole. If, for exam- 
ple, the impacts detected depend on a threshold energy 
of the impinging density wave and the state of the atoms 
making up the detector, then it seems possible that de- 
tection is actually a very rare and, due to the excitations 
of the detector atoms, essentially stochastic event and 
that most density waves remain undetected, but do con- 
tribute to the interference pattern. That such a scenario 



is not as inconceivable as it might seem, can be inferred 
from scanning tunneling microscopy experiments, where 
only a tiny fraction, or 10~ 5 of the actual current in the 
tunnel junction is detected in the experiments [52] . 

In a wider context, there is a long tradition in experi- 
mental physics, starting with Alain Aspect's experiments 
in the 1980s [53], on the non- locality of quantum mechan- 
ics, formalized in the Bell inequalities [54, 55]. These ex- 
periments remain beyond the scope of this presentation, 
although it cannot at present be excluded that non-local 
correlations could play a role in the properties of the 
bivector potential in many-electron systems. Given the 
role of non-local correlations, for example in the con- 
struction of kinetic energy density functionals [56] , such 
a possibility should certainly not be excluded. From a 
different point of view, it can be said that the Bell in- 
equalities [54] arc violated, whenever a system violates 
at least one of following two conditions: (i) Its interac- 
tions are strictly local, and(ii) its spin is counterfactual 
definite (see Rcf. [57] and references therein). Counter- 
factual definite means that an electron possesses a par- 
ticular (= definite) spin, regardless of whether it is mea- 
sured or not. Quantum mechanics, for example, violates 
both of these conditions: the first, in the collapse of the 
wavefunction (which is non-local), the second , in the 
dependence of spin on a measurement (as spin is gen- 
erally isotropic without a measurement). This question 
requires a very detailed and careful analysis, which shall 
be pursued in a follow-up paper. 

From the viewpoint of condensed matter theory and 
computational methods the extension of the equation de- 
rived by Levy, Perdew, and Salmi [9], via a distribution 
of density and field properties, given in Eq. (Ill), is cer- 
tainly the most important aspect of the model. It seems 
at first view surprising that one may describe systems 
composed of a large number of electrons by a wavefunc- 
tion with only four variables, i. e., the amplitudes of den- 
sity and field components, and the direction of the field. 
In particular, since the standard model in many-body 
physics requires at least 3N variables. While the final 
proof that the model is suitable also for many-body sys- 
tems will have to come from the application of the model, 
its predictions, and its agreement with experimental mea- 
surements, it can be analyzed, where the actual gain in 
efficiency could potentially come from. From the view- 
point of DFT, the key problem in orbital free models al- 
ways was the kinetic energy functional [47] . In the Kohn- 
Sham formulation of DFT, this problem is transferred to 
the exchange-correlation functionals, which describe the 
difference in kinetic energy between interacting and non- 
interacting electrons [8]. In the present model, the kinetic 
energy density is given by the Laplace operator acting on 
the density and field components of electrons, e. g. in 
a solid. It is thus different from the von Weizsacker or 
the Thomas-Fermi kinetic energy densities [5]. However, 
as shown in the derivation of the free electron and the 
hydrogen atom, it reverts to these limiting cases with a 
vanishing bivector potential in case of single electrons. 
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In many-body theory, the many-body aspect is usually 
encoded in the construction from single-electron states; 
here, the many-body aspect is accounted for by the bivec- 
tor potential and by the phase coherence of densities and 
field components. The present model of a many-body 
system in a sense does not know of single, i.e., isolated 
electrons, since every single electron is connected to adja- 
cent electrons via coherent density functions and fields. 
In the present model, the many-body aspect of a sys- 
tem is thus encoded in its coherence, and in its com- 
mon energy or chemical potential /x. Since it does not 
know of single electron states, it also does not require to 
build up the many-body wavefunction from Slater deter- 
minants defined via a 3-/V-variable space. It is thus much 
closer to DFT than existing many-body formulations, 
even though electron correlations within solids should be 
accurately described. Considering, that phase coherence 
of also plays a role across the boundary of a peri- 
odic system, most principles known to condensed matter 
theorists which come from periodic boundary conditions 
should be applicable. However, since the potential effi- 
ciency is much higher, in particular in the description of 
large systems, it should also be possible, without lack of 
precision, to simulate very large systems, i.e. mesoscopic 
systems and mesoscopic timescales. 

DFT continues to be a project in development. In 
essence, its successful application to most problems in 
solid state physics and physical chemistry has made it 
an indispensable tool of theorists today. However, the 
current efforts to determine the kinetic energy of a many 
electron system [56, 59] and to go beyond Kohn-Sham 
theory also reveal that the present theoretical framework 
is not an endpoint in the development. One of the mo- 
tivations to develop this theoretical model of electrons 
was the attempt to understand, from the very princi- 
ples of modern physics, why an electron appears such a 
strange entity. The hope was that a better understand- 
ing of fundamental processes and effects may lead to a 
better formulation of the many-electron problem. The 
bivector potential n, introduced in this paper, is a gen- 
uinely novel concept. It arises naturally, as seen, if one 
tries to generalize the LPS equation for a spin system. 
What the exact form of this potential is and how it can 
be formalized in terms of spin-densities, remains to be 
seen. 

We have briefly mentioned in this paper that a rela- 



tivistic theoretical framework, based on the Dirac equa- 
tion and geometric algebra, has been developed by 
Hestenes and others [24, 25]. However, we have not ex- 
tended the present framework to the relativistic domain. 
A reason for this omission is that a non-relativistic frame- 
work seems to be sufficient for the time being, in par- 
ticular in view of applications in condensed matter. In 
this respect it is also worth reflecting the fact that the 
Schrodinger equation, contrary to the wave equation, is 
actually Galilei-invariant (under a kinematic transforma- 
tion of the wave function, see Ref. [34] p. 516). Which 
means that it remains unchanged for systems at low ve- 
locity, which applies to practically all systems studied 
under laboratory conditions. Under these conditions the 
first postulate of the theory of relativity remains thus 
valid, while the second postulate is also valid, since the 
speed of light is much larger than the velocity of elec- 
trons. Imposing an additional Lorcntz-invariance on the 
systems under consideration in DFT seems thus unnec- 
essary at this point. The only real motivation for de- 
veloping such a Lorentz-invariant model of electrons for 
applications in condensed matter physics would be, if the 
description in this manner could be made substantially 
simpler. At present, this seems not the case. However, 
we shall return to this problem in future publications. 
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